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ABSTRACT 

Context. The abundances of the three main isotopes of oxygen are altered in the course of the CNO-cycle. When the first dredge-up 
mixes the burning products to the surface, the nucleosynthesis processes can be probed by measuring oxygen isotopic ratios. 

Aims. By measuring '®0/'^0 and *®0/'®0 in red giants of known mass we compare the isotope ratios with predictions from stellar 
and galactic evolution modelling. 

Methods. Oxygen isotopic ratios were derived from the A'-band spectra of six red giants. The sample red giants are open cluster 
members with known masses of between 1.8 and 4.5 Mq. The abundance determination employs synthetic spectra calculated with the 
COMARCS code. The effect of uncertainties in the nuclear reaction rates, the mixing length, and of a change in the initial abundance 
of the oxygen isotopes was determined by a set of nucleosynthesis and mixing models using the FUNS code. 

Results. The observed *®0/*^0 ratios are in good agreement with the model results, even if the measured values do not present clear 
evidence of a variation with the stellar mass. The observed '®0/'*0 ratios are clearly lower than the predictions from our reference 
model. Variations in nuclear reaction rates and mixing length parameter both have only a very weak effect on the predicted values. 
The '^C/'^C ratios of the K giants studied implies the absence of extra-mixing in these objects. 

Conclusions. A comparison with galactic chemical evolution models indicates that the *®0/**0 abundance ratio underwent a faster 
decrease than predicted. To explain the observed ratios, the most likely scenario is a higher initial *®0 abundance combined with a 
lower initial *^0 abundance. Comparing the measured ratio with the corresponding value for the ISM points towards an initial 

enhancement of as well. Limitations imposed by the observations prevent this from being a conclusive result. 

Key words. Nuclear reactions, nucleosynthesis, abundances - Stars: abundances - Stars: evolution - Stars: late-type 


1. Introduction 


The ratios of the abundances of the three stable isotopes of oxy¬ 
gen, and are important indica tors of the nucle - 

osynthesis and mixing in the interiors of stars (lDearbornlll992l) . 
[n an H-burning environment '®0(p,7)^^F(/S^)^^0 causes a steep 
increase in the abundance until equilibrium is reach ed by 
the competitive *^0(p,a)'^N process (iLandre et alJll99(ih . Ow¬ 
ing to the steep abundance gradient of resulting from in¬ 
complete CNO-cycle burning, the surface abundance of the '^O 
isotope is highly sensitive to the precise depth of convection an d 
the mixing profile in red giants dBoothrovd & Sackmannlll999l). 
jfOis primarily produced via (e.g. lClavtonl 

I 2 OO 3 I 1 . This reaction takes place in the He-burning shells of mas¬ 
sive stars. The isotope is destroyed in hydrogen-burning stars by 
'^0 (p,q')^^N and destroyed during He burning by '^0(a,'y)^^Ne. 
In low and intermediate mass stars, the ^^O abundance thus re¬ 
flects the initial abundance of this isotope and its destruction rate 
during hydrogen burning. The reaction rates for the processes in¬ 
cluding the various oxygen isotopes are relevant for later steps 
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in the stellar nucle osynthesis, in particu lar the s-process (see the 
discussion in, e.g.. lStraniero et an2014l) . 


The solar values for and are 2700 and 

498 (iLodders et al.ll200^ . respectively. Based on the interplay 
between the production, mixing, and destruction of the oxy¬ 
gen isotopes, stellar evolution models predict a correlation be¬ 
tween oxygen isotopic r atios and stellar mass f or red giants 
after the first dredge-up ( Boothroyd et al.|| 1994l: lEl Eidlll994l: 


IStoesz & Herwigl2(K)3l iKarakas & Lattanziol2014t) . The various 

models agree on a steep decrease in the ratio between 

1 and around 2 Mq followed by a weak increase in this ratio for 
higher masses. At the same time, the value after the first 

dredge-up is expected to show a moderate increase for stars be¬ 
low 2 Mq and a constant value afterwards. There are small quan¬ 
titative differences for the predicted ratio between the 

models. A direct comparison of abundances is hampered by vari¬ 
ations in the starting values chosen by the various models. 


For testing the pr edictions, RGB stars of known mass are 
crucial. lEl EidI (1 19941) attempted to do such a comparison with 
the help of a set of bright field stars. However, as noted by the 
author of that study, masses are rather uncertain for these objects. 
Accordingly, a convincing observational test of the models was 
not possible. 


Article number, page 1 of|9] 



























A&A proofs: manuscript no. paper 


The CO vibration-rotation 4.6 pm fundamental or 2.3 pm 
first overtone bands were used in early determinations of O 
isotopes in red giants. The early determinations were limited 
to bright stars due to instrurn ental constraints (see iHarris et~^ 
1 19881 [Smith & Lamberlj[l9^ and references therein). Contem¬ 
porary instrumentation is much more sensitive, allowing the 
study of weaker targets at high spectral resolution and high S/N. 
It is now possible to observe open-cluster (OC) red giants where 
stellar mass can be determined from the cluster age. Exploring 
the causes for oxygen isotopic ratios that deviate significantly 
from predictions based on the solar value provides insight into 
the cosmic matter cycle. 


2. Sample definition and observations 

For the target select ion, we chose the list of red giants in galac¬ 
tic O C provided by iGilroyl (Il989ll plus two stars in NGC 7789 
from jPrugniel et alJ ( 2011 1. The total sample selected on the ba¬ 
sis of observability (see below) and stellar mass includes seven 
stars. iGilrovI (Il989h provides the masses of the stars on the 
giant branch, and we used these values for our analysis. For 
NGC 7789 there ar e two recent a ge determinations of 1.4 Gyr 
(iKalirai et al.l 120081) a nd 1.8 Gyr (Jacobson et all l2Ql ill . When 
using isochrones from iBressan et dJ (I2012 ). the turn-off mass is 
1.8+0.2 M(7). The metallicities of the select ed clusters are close 
to solar (iGilrovll 19891 : jPrugniel et al.ll20l'll) . Our sample of red 
giants is listed in Table[T] 

Subsequent to taking the observations, we found that 
HD 49050 is no longe r considered a member of NGC 2287 
(iMermilliod et al J2008l) . On the other hand, cluster membership 
of HD 16068 in Tr 2 a nd HD 27292 in NGC 1545 was confirmed 
by Zacs et al.l (1201 ll) . For HD 68879. iFrinchabov & Maiewskil 
(I 2 OO 8 I 1 list a cluster membership probability of more than 96 %. 
For the other stars, their radial velocities are very clos e to the 
respective mean cluster velocity (IMermilliod et al.ll2008[) . which 
supports their membership. A re-examination of the kinematics 
of the UMa group based o n Hipparcos data leaves the member¬ 
ship of HD 30834 unclear (iKing et alj|2003h . 

We obtained several small pieces of high-resolution spec- 
tra in the H and K bands using the Phoenix spectrograph 
dHinkle et alj|200(il) at the Kitt Peak 2.1 m telescope. Observa¬ 
tions were obtained in December 2013 and January 2014. The 
standard infrared observing procedure with two nodding posi¬ 
tions was applied. Telluric lines were removed by ratioing the 
spectra to spectra of hot stars observed at similar airmass. The 
telluric lines in the hot star spectra were also used for the wave¬ 
length calibration. 


3. Isotopic ratios 

Our goal was to derive the ^®0/*^0, '®0/'^0, and '^C/'^C ra¬ 
tios for our seven target stars. Owing to the very limited spectral 
range of about 100 A covered by a single observation, the num¬ 
ber of spectral features available for our analysis is smaller than 
in the earlier studies based on FTS scans of the whole /f-band 
region. To measure the ^^O abundance we used two to four 2- 
0 band C^^O lines located near 4285.4 cm *. These lines are 
largely unblended and can be easily identified for a stellar tem¬ 
perature of up to 4200 K. 

Determining the '^O abundance from JiT-band spectra is more 
difficult. We ultimately used a region around 4226 cm ', which 
was used by iGarcla-Hernandez et al.l (1201 Oil for their study of 
oxygen isotopes in RCrB stars. The most usable line of C'^O in 


this region is the 2-0 R23 line, which is not affected significantly 
by telluric lines. However, this line is affected by an unidentified 
blend, so its usability for abundance determination is limited to 
7’eff<4100K and logjg<2.0. An attemp t to use the 2-0 band head 
ofC'^O, as in ISmith & Famber3(ll990l) . for instance, failed since 
the band head is located too close to the transmission cut-off of 
the Phoenix order separation filter. 

Complimentary information on mixing in the stellar atmo¬ 
sphere was derived by measuring the ratio for our sam¬ 
ple stars. Several first overtone and lines near 

4246 cm“' and 4262 cm“' were used. We also attempted ob¬ 
servations of the CO fundamental spectrum near 4.6pm where 
many CO lines from all three oxygen isotopes can be found. 
However, lower flux levels, line blending, and difficulties in set¬ 
ting the continuum level limited the use of these data. The data 
presented in this paper are part of a larger observing programme 
on oxygen abundances in cool red giants. Details will be pre¬ 
sented elsewhere. 


Isotopic ratios were computed using spectrum synthesis 
techniques. For the stellar atmospheric structure, we used 
the hydrosta t ic sph erical COMARCS models described in 
lAringer et al.l (l2009h . Synthetic spectra were calculated with the 
COMA code using the same set of opacity data. Calculations 
were done under the assumpti on of FTE. We adop ted values for 
solar composition provided bv ICaffau et al.l (l2009h . Abundances 
were determined by a direct comparison between observed and 
synthetic spectra using both visual inspection and the measured 
line depths. For a more detailed description of our approach, we 
refer to earlier applic ations to the fitt i ng of simi lar high reso¬ 
lution spectra given in iFebzelter et al.l (l2008l) and iFederer et al.l 
(l2009ll . Fine positions throughout the st udied spectra l range have 
been improved b y using the work of iHinkle et al.l (1 19951) and 
iHase et akl (I 2 OIOI) . Data for the lines used in this analysis are 
listed in Table lA.il 


The line strengths depend on the stellar temperature, surface 
gravity, chemical composition, and abundance of each isotope. 
For the stars taken from the list of iGilrov ( 19891) .we u s ed the 
values for and logg given there. McDonald et akl (12012l) 
independently determined Teff values for HD 27292 (3844K), 
HD 30834 (4247 K), and HD 68879 (4552 K), which are in rea¬ 
sonable agreement with the values given by GilrojQ. For the 
two stars in NGC 7789, the s tellar parameters were taken from 
the study of iPrugniel et al.Jj 20111). Sin c e all clusters have a 
metallicity close t o solar (lGilrovTl989t Llacobson et alJ 1201 ll 
iMilone et al.ll201 ll) . we calculated our grid of synthetic spectra 
only for solar composition. For C/O we chose a fixed value of 
0.3 to resemble a typical value for post first dredge-up composi¬ 
tion. The microturbulence ^ was s et 2.0 km s~*, w hich is a typical 
value for these stars according to lGilrovI (Il989l) . The macro tur¬ 
bulence was set to values between 1.5 and 2.5 kms^' to optimize 
the fit of all stellar lines within the observed wavelength range. 
All parameters taken from the literature were cross-checked and 
confirmed with our spectra. The finally chosen values for Teff 
and log g are given in columns 4 and 5 of Table[T] Examples for 
the spectral fit are shown in Figs. [1] and |2 We note that our syn¬ 
thetic spectra cannot fit the cores of the strongest CO lines (low¬ 
est excitation) properly (e.g. Figl B. This diffic ulty has also been 
encountered by other authors (e.g. lTsuiill2008h and likely results 


* ICesetti efd] (120131) list a rather high temperatur e of 6427 K for 
HD 16 068. They refer to a parameter determination bv iHolmberg et all 
( I 2 OO 8 I) . However, we could not find this star in the la t ter ca talogue. 
Therefore, we suppose there is a mistake in ICesetti et all (1201.31) . 
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Fig. 1. Observed spectrum of HD 27292 (dots) showing two '^C'^O 
lines. The three synthetic spectra (solid lines) are for 7’eff=4000K and 
logg=1.5 with ‘^0/‘’0=100, 300, and 500. 

from limitations of hydrostatic model atmospheres in describing 
the outer layers of a red giant properly. 

The resulting oxygen isotopic ratios are listed in Table[T] The 
two stars in the cluster NGC7789 show similar values. For the 
hottest star in our sample, HD 68879, we could not determine 
and because the CO lines were too weak. For 

HD 30834 we could not measure because the contri¬ 

bution of the C'^O line to the blend could not be constrained 
properly. For this star, the value given for in Table[T]is 

a lower limit. We also determined carbon isotopic ratios for all 
our stars. The carbon isot ope results are in good agreement with 
the values determined bv lGilro\^(l 19891) (see Table[T]i. 

To interpret the and ratios, we need to check 

whether a solar abundance value is appropriate for our sam¬ 
ple stars. We did this by computing model spectra with altered 
oxygen abundances and with stellar parameters and carbon iso¬ 
topic ratio as derived in the previous step. Since C/O < 1 the 
strengths of the CO lines are dependent on the carbon abun¬ 
dance and the isotopic ratios of oxygen but not on the total oxy¬ 
gen abundance. To derive the oxygen abundance, OH vibration- 
rotation lines in the H band were modelled. The //-band spec¬ 
tra around 6072 cm ' were obtained for six of the K giants 
in Table[T] We modelled the three least blended '®OH lines, 
4-2Pif6.5, 2-0Pifl6.5, and 2-0P 2 e 15.5. The resulting abun¬ 
dances listed in Tab le[T] were all sub-solar by -0.16 to -0.36 dex. 
The uncertainties are dominated by the temperature sensitivity 
of the chosen OH lines. As noted above, the total oxygen abun¬ 
dance has no effect on the strengths of the CO lines. We also 
do not expect that the oxygen abundance will affect the atmo¬ 
spheric structure because photospheric H 2 O is not present in the 
temperature rang e of our sampl e stars. 

According to iGilrovI (Il989l) . the uncertainties of the stellar 
parameters are A Teff +150K, Alogg +0.3, +0.2kms“', and 

A [Fe/H] +0.2. Alternative parameter values found in the litera¬ 
ture agree within these error bars. For each star we derived the 
isotopic ratios by changing the stellar parameters within these 
ranges. Some combinations, however, were not included when 
estimating the uncertainties because the corresponding synthetic 
spectra do not provide a good fit to the observations. Since the 
continuum level in the K band is well defined for these stars, we 



2.3655 2.3660 2.3665 

X [a'^] 


Fig. 2. Observed spectrum of HD 27292 (dots) with a line 

and a line marked. The three synthetic spectra (solid lines) are 

for reff=4000K and logg=1.5 with ‘®O/‘*O=300, 500, and 700 and 
*^C/*'^C=15 and 30, respectively. 

did not consider it to have an impact on the error budget. The 
final uncertainties, given in Table[T] were then determined by 
combining the maximum differences from changes in the stel¬ 
lar parameters with the scatter resulting when various lines were 
used for the abundance determination. The latter was typically 
a factor of 2 smaller than the uncertainties from the stellar pa¬ 
rameters. As noted above we did not determine the C/O ratios 
from our spectra but set them to a fixed value of 0.3 in our anal¬ 
ysis. We explored the effect of changing the C/O ratio by +0.1 
dex on the derived oxygen isotopic ratios. Such a change some¬ 
what modifies the strength of all CO lines, but primarily those 
of weak or moderate strength. To achieve a good fit of the whole 
spectrum again, the stellar parameters have to be adapted so that 
the net effect on the isotopic ratios is small and covered by the 
error budget derived above. The effect of changing the macro¬ 
turbulence value by 1 kms“' is a few percent on the resulting 
isotopic ratio. 

For HD 16068 and HD 27292, we also observed the M-band 
spectrum near 2154cm“*. As mentioned above, the M-band 
spectra could not be used to derive reliable isotopic ratios. How¬ 
ever, attempting to achieve a reasonable fit to the 4.6 /rm spectra 
limits the ^®0/'^0 and '®0/'®0 ratios to a range between 200 and 
600, in agreement with the findings from the K band. 

4. Nucleosynthesis and mixing models 

To interpret the measured oxygen isotopic ratios, we have com¬ 
puted a set of stellar models by means of the FUN S (FUll Net¬ 
work Stellar evolution) code (IStraniero et al.ll200^ . In practice, 
all the stable isotopes from H to F have been explicitly included 
in the H-burning nuclear network coupled to the stel lar structure 
equati ons. The reaction rates have been taken from llliadis et al.l 
mM and lAdelberger et ^ (l201 ll) . The mixing length parame¬ 
ter, a'=1.82, was calibrated by means of a standard solar model 
computed with the same set of nuclear reaction rates, equa¬ 
tion of state, radiative opacities, and composition as the models 
used in this study. Initial metallicity (Z=0.015) and He abun¬ 
dance (Y=0.27) are the early solar values as derived from the 
same standard solar model. The solar composition was taken 
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Table 1. Stellar parameters and isotopic ratios. Data in columns 3 to 6 are taken from the literature. 


Cluster 

Star 

Mass [Mq] 

Teff [K] 

log^ 

‘^C/‘^C Lit. 

Ref. 

[‘'’O/H] 

Ibo/l/Q 

Ibo/ISQ 

^^crc 

NGC 7789 

971 

1.8+0.2 

3700 

1.2 

- 

(2) 

-0.18+0.1 

300!^“ 

340!^^° 

15+8 


415 

1.8+0.2 

3800 

1.2 

- 

(2) 

-0.16+0.1 

2501^“ 

270+150 

19+8 

NGC 2548 

HD 68879 

2.6+0.3 

4600 

2.2 

25 

(1) 

- 

- 

- 

16+7 

NGC 1545 

HD 27292 

2.8+0.2 

4000 

1.5 

25 

(1) 

-0.33+0.1 

240!^° 

340!*^° 

27+10 

UMa Grp. 

HD 30834 

2.9+0.4 

4000 

2.2 

26 

(1) 

-0.36+0.1 

260^^^ 

z.uu_50 

>200 

22+7 

NGC 2287 

HD 49050 

_a 

4100 

1.6 

23 

(1) 

-0.19+0.1 

300!^° 

J} 

30: 

Tr2 

HD 16068 

4.5+0.4 

4000 

1.5 

30 

(1) 

-0.28+0.1 

280+®° 

370+100 
j /u po 

28+10 


References. fll lGilrot^ ( 119891) : 12) lPrugniel et al.l (1201 Ih . 

Notes. Star is probably not a cluster member. S/N at the wavelength of the line too low. 


from th e list provided by Lodders et al.l (|2009|) . For details we 
refer to iPiersanti et al.l (12007 ). A comparison of FUNS predic- 
tio ns on stellar m asses and cluster ages with the values derived 
by iGilrovI (Il989h showed good agreement within the uncertain¬ 
ties, so we refrained from re-determining the stellar masses of 
our sample stars. We point out that there is a minor inconsis- 
ten cy in our study in the sense that the oxygen abundance given 
by iLodders et al. ( 20091) is 0.07 dex lower than the value of 
ICaffau et ah (2009) we used for computing the model spectra. 
Since the difference is very small, this inconsistency has negligi¬ 
ble consequences on our analysis. 


Stellar masses between 1.8 and 5.0Mq have been investi¬ 
gated. The complete set of model parameter combinations and 
the corresponding results are summarized in Table|3 The first 
block (first 6 rows) refers to the reference models (R), those 
obtained with the recommended values of the reaction rates, 
the calibrated mixing length, and the solar composition. Then, 
to quantify the theoretical uncertainties, additional models have 
been computed by varying nucleosynthesis and mixing inputs, 
as well as the initial abundances of the three stable O isotopes. 


As usual, all the models presented here include a treatment 
of the convective mixing. Therefore, red giant models show the 
composition modified by the first dredge-up. No extra mixing in¬ 
duced, for example, by rotation, thermohaline circulation, mag¬ 
netic buoyancy, or gravity wave has been considered. As is well 
known these processes are hampered in red giant stars by the 
sharp molecular weight gradient left by the first dredge-up. Only 
in stars with M < 2 Mq does the shell H-burning reach the H 
discontinuity during the RGB phase, thus smoothing down the 
p - gradient. With the possible exception of the two NGC 7789 
stars, all the other stars in our sample have masses higher than 2 
Mq. In the model with M= 1.8 Mq, the H-burning attains the H 
discontinuity when log(L/Lo) ~ 2. The two giant stars observed 
in NGC 7789 are slightly brighter than this thresholcfl so some 
extramixing has possibly modified their compositions. However, 
the observed (Table[T]) is only slightly smaller than the 

value expected after the FDU and, in any case, within the error 
bar. The O is otopic ratios are n ot affected by moderate extramix¬ 
ing (see, e.g.. lAbia et al.ll201^ . 


^ The luminosity may be estimated by means of Tcff and log g listed 
in Tabled 


5. Discussion 

5 . 1 . 1^0 

All of the six sample K-giants with determined oxygen isotopic 
ratios show a very similar ratio between 250 and 300. 

Reference models, i.e. those obtained starting from a solar com¬ 
position, account for the observed values within the observa¬ 
tional errors. On the other hand, this agreement does not nec¬ 
essarily imply that the protostellar abundances were nearly 
solar. In the H-burning shell, the is directly linked to the 
through the NO cycle (see Figure|3]l. Therefore, the final '®0/*^0 
is basically fixed by the initial abundance of the most abundant 
isotope of the cycle, '®0, with variations of the initial '^O having 
a negligible effect. This is demonstrated by model C170H with 
an initial ’^O abundance 34% higher than the reference model 
but nonetheless a final '®0/*^0 ratio that is practically identical 
to the reference model. However, a variation in the initial '®0 
will modify the isotope ratios after the first dredge-up. Never¬ 
theless, a partial redistribution of the '®0 excess or deficiency 
over all the isotopes of the NO cycle mitigates the change in the 
final '®0/'’0. All the six giants in our sample show clear sig¬ 
natures of subsolar ['®0/H], with an average value -0.25 + 0.03 
(see TablellJ. Models C160L demonstrate that such a 40% re¬ 
duction of the initial would imply a ~ 14% reduction of the 
final ratios, which is still compatible with the observed 

ratios (within Icr; see also the black-solid curve in Figure|4]). 

Variations in the nuclear reaction rates have weak effects 
compared to the observational uncertainties. A change of Z from 
0.15 to 0.02 implies a 10% increase in the '®0/'’0. On the con¬ 
trary, a 10% reduction is obtained if Y is changed from 0.27 to 
0.32. 

As a whole, at variance with model predictions, the measured 
ratios do not show any clear evidence of a variation with 
the stellar mass. However, the small number of observations, also 
affected by a rather large error, hampers any more thorough in - 
vestigation of this issue. As outlined in ILebzelter et al.l (12012h . 
reducing the errors is hindered by fundamental problems in de¬ 
termining stellar parameters for cool giants. Effects of changes 
in the mixing length parameter are similarly far below the obser¬ 
vational detectability. 

5 . 2 . 1*^0 

The models show that the ratio after the first dredge-up 

is almost constant (for M> 2Mo). However, the five K-giants 
analysed here have typical values near 350, which is about half 
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Table 2. Effects of theoretical uncertainties on and *®0/**0. 


Initial isot. ratios Post-FDU isot. ratios 


label 

M 

(Mo) 

160/170 

160/180 

'^O-i-p 

(rates) 

'^o+p 

(rates) 

^ml 

12c/13cfl 

160/170 

160/180 

Reference models^’ 

R 

1.8 

2696 

499 

recom. 

recom. 

1.82 

24.5 

440 

684 

R 

2.0 

2696 

499 

recom. 

recom. 

1.82 

24.2 

262 

703 

R 

2.5 

2696 

499 

recom. 

recom. 

1.82 

23.8 

208 

711 

R 

3.0 

2696 

499 

recom. 

recom. 

1.82 

23.5 

253 

709 

R 

4.0 

2696 

499 

recom. 

recom. 

1.82 

23.2 

338 

703 

R 

5.0 

2696 

499 

recom. 

recom. 

1.82 

22.9 

381 

704 


Nuclear reactions changed^ 


017L 

4.0 

2696 

499 

low 

recom. 

1.82 

23.2 

320 

704 

017H 

4.0 

2696 

499 

high 

recom. 

1.82 

23.2 

355 

703 

017L 

2.5 

2696 

499 

low 

recom. 

1.82 

23.8 

193 

711 

017H 

2.5 

2696 

499 

high 

recom. 

1.82 

23.8 

224 

711 

018L 

4.0 

2696 

499 

recom. 

low 

1.82 

23.2 

336 

681 

018H 

4.0 

2696 

499 

recom. 

high 

1.82 

23.2 

337 

722 

018L 

2.5 

2696 

499 

recom. 

low 

1.82 

23.8 

210 

693 

018H 

2.5 

2696 

499 

recom. 

high 

1.82 

23.8 

210 

727 

BL 

2.5 

2696 

499 

hi( 7 ),lo(a) 

low 

1.82 

23.8 

198 

693 


Initial composition changed"' 


C170H 

2.5 

2012 

499 

recom. 

recom. 

1.82 

23.8 

203 

711 

C180H 

2.5 

2696 

332 

recom. 

recom. 

1.82 

23.8 

208 

474 

C180HH 

2.5 

2696 

249 

recom. 

recom. 

1.82 

23.8 

209 

356 

C180HH 

3.0 

2696 

249 

recom. 

recom. 

1.82 

23.5 

253 

355 

C180HH 

4.0 

2696 

249 

recom. 

recom. 

1.82 

23.2 

336 

352 

C160LL 

2.5 

1348 

249 

recom. 

recom. 

1 82 

23.7 

169 

362 

C160LL 

3.0 

1348 

249 

recom. 

recom. 

1.82 

23.5 

203 

361 

C160LL 

4.0 

1348 

249 

recom. 

recom. 

1.82 

23.1 

258 

359 

C160L 

1.8 

1617 

299 

recom. 

recom. 

1.82 

24.5 

343 

413 

C160L 

2.0 

1617 

299 

recom. 

recom. 

1.82 

24.1 

216 

426 

C160L 

2.5 

1617 

299 

recom. 

recom. 

1.82 

23.8 

180 

432 

C160L 

3.0 

1617 

299 

recom. 

recom. 

1.82 

23.5 

217 

431 

C160L 

4.0 

1617 

299 

recom. 

recom. 

1.82 

23.1 

278 

428 

C160L 

5.0 

1617 

299 

recom. 

recom. 

1.82 

22.9 

307 

429 


Mixing length parameters a^i changed 


AL 

AH 

2.5 

2.5 

2696 

2696 

499 

499 

recom. 

recom. 

recom. 1.50 
recom. 2.00 

23.8 

23.8 

208 

208 

711 

711 

Metallicity change: Z=0.02 

Zvar 

2.5 

2696 

499 

recom. 

recom. 1.82 

23.6 

189 

723 

He change: Y=0.32 

Yvar 

2.5 

2696 

499 

recom. 

recom. 1.82 

23.8 

232 

709 


Notes. Initial alwa ys set to 89. Reference models as obtained by assuming solar composition dLodders et alJl20Q^ . 0-t-p reaction 

rates from llliadis et_alj j20K]h and so lar-calibrated mixing length. Models obtained using the lower and upper proton capture rates of or 
suggested bv llliadis et all (l2010h . All other parameters equal to reference models. Model BL: capture rates setting chosen to minimize **0 
depletion. ^ Initial composition changed: ‘^0x1.34 (C170H), '®Oxl.5 (C180H), ‘*0x2 (C180HH), ‘«Ox0.5 (C160LL), ‘«Ox0.4 (C160L). 


the value predicted by the reference models (R). None of the 
sample stars has an abundance in the literature. While our 
values seem to be amongst the lowest ratios measured 

in evolved giants, values around 400 have been mea¬ 

sured in field giants. [Harris et’aP (Il988h report similar values 
for the bright stars a Ari, aSer, jSAnd, ySPeg, and /lUMi. The 
[Harris et"^ (Il988t) study used different molecular lines than we 
have selected, which indicate that such low ratios are probably 
not the result of an inapprop riate selection of sp ectral lines. All 
but one red giant studied by [Harris et"^ (Il988h were found to 
have ratios <600. 


As for ‘^O, our tests of the uncertainties of the nucleosyn¬ 
thesis model show that the only route to significantly modifying 
the ratio after the first dredge-up is to change the initial 

abundances of individual isotopes. At variance with the '^O, for 
which only a variation in the initial abundance affects the 
final a modification of the initial could play an im¬ 
portant role in interpreting the measured ratio. Owing 

to the weakness of the '^0(/7 ,reaction, remains de¬ 
coupled from the CNO cycle (see Figure^. On the other hand, 
it can be easily destroyed by fast proton captures (mainly by 
the ‘^0 (p,q')^^N reaction) so that the '®0/^®0 ratio should in- 
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Fig. 3. Illustration of the network of nuclear reactions within the CNO 
cycle. 


crease after the first dredge-ujfl Our models show that for any 
choice of the initial abundance pair and for any initial 

mass, the final will always be ~ 1.4 times greater than 

the corresponding initial value. For this reason, starting from a 
solar value of the ratio, ~ 500, it is impossible to ac¬ 

count for the measured values. Our models require that either the 
'®0 or the '*0 be varied in order to match the observations. For 
instance, model C180HH with an initial '^O twice solar abun¬ 
dance matches the observed '®0/'*0. However, this is also the 
result when using model Cl60LL with an initial '®0 abundance 
half solar. It follows that combinations of these two assumptions 
also provide a good reproduction of the observations. For in¬ 
stance, assuming a 40% reduction of the initial '®0, as implied 
by the observed average ['®0/H], results in which is in 

reasonable agreement with the measured value (within Icr; Fig- 
ure|5]). The agreement can be improved further by increasing the 
initial about 10%. 

5.3. Implications of the results 

By combining ['®0/H], '®0/’^0, and measurements in 

six giants belonging to galactic OC of known ages, we can infer 
the O isotopic composition of the respective protostellar nebu¬ 
lae. In particular, the optimal protostellar abundances modelled 
from the observed stellar abundances imply subsolar '®0, i.e. 
-40% on the average, and slightly supersolar i.e. H-10%. 
The ratios in the protostellar gas are reduced by a fac¬ 

tor of ~ 1.4 from the values measured in the red giant stars. No 
direct information on the protostellar can be obtained from 
our observations. 

Is this picture that arises from our measurements in OC gi¬ 
ants compatible with the predictions of extant models of galac¬ 
tic chemical evolutiorQ (GCE)? is a primary product of the 
He burning. It is produced through the main He-burning chain; 
3a —>^^C-i-a is also produced in He-burning regions 

as a secondary product synthesized through the following chain: 
^"^N(a, 7 )^^F(j 6 )^® 0 , where the abundance of is equivalent, 
in practice, to the original amount of C-i-N-i-O. The most effi¬ 
cient polluters of both and are massive stars exploding 
as core-collapse supernovae. is a secondary product of the 

^ As previously noted is practically unaffected by the first dredge- 
up. 

Isotopic ratios of oxygen measured in stars are typically given as 
liQjxxQ jjjg literature. Accordingly, we used this format of the iso¬ 
topic ratio in the first part of the paper. Galactic chemical evolution 
models, however, often give the reversed ratio '“0/*®0. Therefore, we 
decided to use “O/'^O in this section to allow for an easier comparison 
with the literature. 



Fig. 4. ratios after the first dredge-up versus stellar mass in Mq. 

Lines represent the theoretical predictions, and the red squares are ob¬ 
servations. Blue solid line: reference model (R). Blue dashed lines: 
proton capture rate modified within the suggested upper and lower rates 
(017L and 017H). Solid black line: theoretical predictions as obtained 
by reducing the initial ‘'^O to [‘'^O/H]=-0.22 (C160L). 


H burning (CNO cycle, see Figure O. Therefore, in addition to 
massive stars, the winds of intermediate and low mass stars con¬ 
tribute on a longer timescale to the oxygen pollution of the ISM 
by contributing '^O. In this framework, it is not a surprise that 
GCE models generally predict an increase with time in both the 
160 /I 80 and the as expected for the ratio of a primary 

to a secondary product. Less straightforward is the evolution of 
the ratio of the two secondary O isotopes. In general, when two 
secondary isotopes are produced by the same class of stars, their 
ratio should remain constant. However, this is not the case for 
which is expected to decrease when low and intermedi¬ 
ate mass stars begin to contribute to the pollution of the lightest 


ate mass stars bepn to contribute to tne pollution ot me ligntest 
isotope (C lavtonlTl98i^ iPrantzos et alJ 19961: iMever et al.ll2008t 


iKobavashi et al.ll201 itlNittler & Ga d^ 20H 

In Eigure |6] predictions of the Prantzos et alJ (1 19961) GCE 
model are compared with our estimates of in the parent 

nebulae of the OC. The solid line represents the predicted evo¬ 
lution of the ratio over the past 6 Gyr. The CGE model 

refers to abundances in the interstellar medium (ISM) at the cur¬ 
rent solar galactocentric radius, Rqc - 8.5 kpc. The ra¬ 

tio is plotted relative to the solar value while time is scaled to 
the epoch of the solar system formation. Red squares show our 
isotope ratio estimates for the three clusters NGC 7789, NGC 
1545, and Tr 2. In spite of the large uncertainties, it appears 
that the primordial gas of the three clusters underwent a faster 
increase in the abundance than predicted from the GCE 
model. This figure do es not change by using the more recent 
IKobavashi et al.l d201 ih GCE model. The galactocentric distance 
of the three clusters is only slightly greater than solar, between 
9kpc and 9.7kpc. Although migration of the Sun and the clus- 
ter0 from small galactocentric radius cannot be excluded, the 
birthplaces should be within 6 < Rcc < 10kpc. 

Oxygen isotopic ratios in molecular clouds located in this 
portion of the Galaxy are available from infrared and ra- 


^ Effects of migration are probably small for the two younger clusters, 
NGC 1545 and Tr 2, while it cannot be excluded in the case of NGC 
7789 (1.7 Gyr old). 
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Fig. 5. ratios after the first dredge-up versus stellar mass in 

Mq. Lines represent the theoretical predictions, and the red squares are 
observations. Solid blue line: reference model (R). Dashed blue lines: 
'*0 proton capture rate modified within the suggested upper and lower 
rates (018L and 018H). Solid black line: model predictions obtained 
by reducing the initial '«0 to [‘^O/H]=-0.22 (C160L). Dashed black 
line: model with the initial ’*0 abundance increased by a factor of 2 
(C180HH). 


dio observations ( Wouterloot et al. 20081 Milam et’aP 120051 : 
iPolehampton et all 2005 1. Nittler & GaidosI ( 2012h ha\^derived 
iSQ/ifio for molecular clouds by combining re¬ 
ported bv lWouterloot et al] (12008 ) with the galactic gra¬ 

dient obtained bv iMilam et al.l ( 2005h from CO observations. 
The two triangles in Figure |6] represent the weighted averages 
of these measurements, as obtained from two different molecu¬ 
lar lines, CO J=l-0 and CO J=2-l. Only molecular clouds for 
6 < Rgc < lOkpc hav e been considered. Other measurements, 
derived fro m OH lines (^lehampton et alJl2005h and formalde¬ 
hyde data dWilson & Roodlll994 ), are within the range covered 
by the CO measurements. In the same figure, the shaded area 
is representative of the whole set of available measurements for 
the ISM lying on the galacti c disk approximately be tween 6 and 
lOkpc (see also figure 2 in iNittler & GaidosI l2012h . Compared 
to the ISM, the '®0/^®0 ratios we have derived for three OCs are 
close to the upper bound of the available measurements. 

Although we cannot constrain the protostellar '^O from 
OC giant observations, molecular clouds located between 6 
and 1 kpc from the galactic centre show that the ra¬ 

tio ra nges between 3 a nd 5, with an a verage v a lue 4 .16 + 
0.09 (IW outerloot et al.l (l2008h : see also iPenziasI (Il981h and 
iMilam et alJ ( 2005 )1. Since the ratio of these two secondary iso¬ 
topes can either remain constant or decrease with time and since 
our measurements indicate slightly supersolar in the ISM 
from which the host clusters were born, a similar or even larger 
enhancement of is expected, otherwise the initial 
would be definitely too large compared to the typical values 
found in the galactic disk. 


6. Conclusions 

We measured the ratios of the three main isotopes of oxygen in 
the atmospheres of six evolved K giants. Except for one the stars 
are members of stellar clusters and therefore have well-defined 



time (Gyr) 

Fig. 6 . Evolution of **0/*®0 with time. The isotopic ratio is plotted 
relative to the solar value. The zero point of the time axis is defined by 
the formation of the sun. The solid black line indicates the predictions 
from galactic chemical evolution (GCE) models. Red boxes mark our 
estimates of in the parent nebulae of the three OCs NGC 7789, 

NGC 1545, and Tr2. The two measurements for NGC 7789 have been 
combined to one point. Blue symbols indicate ISM values. The shaded 
area is representative of the whole range of ISM measurements. 

masses. The and ratios are compared with 

model predictions. There is agreement between observations and 
predictions for the ratio. The dependency of on 

mass, expected from nucleosynthesis models, is not observed. 
However, the predicted variation is comparable to the data er¬ 
rors. The observed ratio can be brought into agreement 

with model predictions if the initial isotopic ratio is about half of 
the solar value. Uncertainties in proton capture rates and mixing 
length cannot account for the observed difference. By combin¬ 
ing this result with our estimation of the ['®0/H], we conclude 
that a moderate enhancement of the initial abundance rela¬ 
tive to the solar value in combination with a subsolar abun¬ 
dance can provide a good match to the observations. Finally, to 
maintain the protostellar in agreement with the values 

typically measured in nearby molecular clouds, we also infer a 
moderate enhancement of the initial abundance. 
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Appendix A: Line data 
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Table A.l. Molecular lines of (selection), *^C**0, and *®OH used in the analysis. 


Molecule" 

Transition 

wavelength* 

[A] 

E* 

[cm“'] 

gf 

12cl6o 

2-0 R1 

23432.693 

3.845 

1.802E-07 


4-2 R41 

23540.840 

7492.257 

2.962E-05 


4-2 R60 

23547.324 

11086.589 

4.871E-05 


4-2 R38 

23553.885 

7043.948 

2.698E-05 


4-2 R24 

23664.988 

5390.375 

1.583E-05 


3-1 R4 

23667.714 

2181.369 

1.392E-06 


2-0 P9 

23680.105 

172.978 

7.605E-07 


2-0 R14 

23502.671 

393.276 

1.407E-06 


2-0 R30 

23316.224 

1737.760 

3.223E-06 


2-0 R29 

23324.897 

1625.950 

3.099E-06 


2-0 R28 

23333.966 

1517.827 

2.976E-06 


2-0 R26 

23353.313 

1312.658 

2.735E-06 

12cl8o 

2-0 R23 

23656.947 

1009.011 

2.298E-06 

13c16o 

2-0 R68 

23501.045 

8501.143 

8.980E-06 


2-0 R69 

23507.818 

8747.480 

9.172E-06 


2-0 R70 

23515.015 

8997.172 

9.367E-06 


2-0 R71 

23522.680 

9250.210 

9.565E-06 


2-0 R76 

23567.025 

10565.262 

1.059E-05 


2-0 R19 

23660.695 

697.622 

1.878E-06 

'®OH 

4-2Piy6.5 

16477.325 

7683.8196 

1.800E-05 


2-0 Pi/ 16.5 

16460.532 

4915.1296 

1.114E-05 


2-0 Pi,, 15.5 

16452.548 

4939.7718 

1.047E-05 


Notes. CO line data from lGoorvitchl dl994h : OH lines from iRothman et aO ( 120091) ® vacuum lower state term energy 
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